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Abstract We have collected all available spectra and photometric data from SDSS cat- 
alogue for bright AGNs complied from the first three months of the Fermi large area 
telescope all-sky survey. Based on the 106 high-confidence and 11 low-confidence asso- 
ciation bright AGN list, the photometry data are collected from SDSS DR7 for 28 sources 
(12 BL Lacs and 16 FSRQs), two of which are low-confidence association bright AGNs. 
Among these 28 SDSS photometric sources, SDSS spectra are available for 20 sources (6 
BL Lacs and 14 FSRQs). The black hole mass Mbh and the broad line region (BLR) lu- 
minosity were obtained for 14 FSRQs by measuring the line-width and strength of broad 
emission lines from SDSS spectra. The broad emission lines measurements of five FSRQs 
are presented for the first time in this work. The optical continuum emission of these 14 
FSRQs is found to be likely dominated by the nonthermal jet emission through comparing 
the relationship between the broad Mg II line and continuum luminosity to that of radio 
quiet AGNs. The black hole mass of 14 FSRQs ranges from to lO^^M©, with 

most of sources larger than 1O^M0. The Eddington ratio iboi/iEdd ranges from 10~^'^ 
to ^ 1. This implies that the optically thin, geometrically thick accretion disk may exist 
in these FSRQs. 

Key words: galaxies: active — galaxies: nuclei — gamma rays: observations — quasars: 
emission lines — BL Lacertae objects: general 



1 INTRODUCTION 

Blazars, including BL Lac objects and flat-spectrum radio quasars (FSRQs), are the most extreme class 
of active galactic nuclei (AGNs), characterized by strong and rapid variability, high polarization, and 
apparent superluminal motion. These extreme properties are generally interpreted as a consequence of 
non-thermal emission from a relativistic jet oriented close to the line of sight. As such, they represent 
a fortuitous natural laboratory with which to study the physical properties of jets, and, ultimately, the 
mechanisms of energy extraction from the central supermassive black holes. 

The most prominent characteristic of the overall spectral energy distribution (SED) of blazars is 
the double-peak structure with two broad spectral components. The first, lower frequency component 
is generally interpreted as being due to synchrotron emission, and the second, higher frequency one as 
being due to inverse Compton emission. BL Lac objects (BL Lacs) usually have no or only very weak 
emission lines, but have a strong highly variable and polarized non-thermal continuum emission ranging 
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from radio to 7-ray band, and their jets have synchrotron peak frequencies ranging from IR/optical to 
UV/soft-X-ray energies. Compared to BL Lacs, FSRQs have strong narrow and broad emission Hnes, 
however generally have low synchrotron peak frequency. According to the synchrotron peak frequency, 
BL Lac objects can be divided into three subclasses, i.e. low frequency peaked BL Lac objects (LBL), 
intermediate objects (IBL) and high frequency peaked BL Lac objects (HBL) (Padovani 8c Giommi 
1995). 

Although a point of agreement about the second component is that the 7-rays of blazars are produced 
in relativistic jets by inverse Compton scattering, the origin of the seed photons (optical/IR), the location 
and size of the emitting region, and the degree of relativistic beaming of the high-energy radiation, are 
all still unknown. Seed photons may come from the jet itself (Synchrotron Self-Compton model: SSC), 
from an accretion disk around a supermassive black hole at the base of the jet, or else from photons of 
the broad emission line region (e.g. Wehrle et al. 1998; Collmar et al. 2000). If the 7-rays origins from 
up-scattering of broad emission line photons, a correlation between the 7-rays and the emission lines 
would be expected. However, a direct comparison shown that there is no evidence of correlation between 
the 7-rays and broad emission lines luminosity (Fan 2000), likely due to the lack of line data. Through 
estimating the energy density in the relativistic jet blobs contributed by the broad line region (BLR), Fan, 
Cao & Gu (2006) claimed that the 7-ray emission can be from external Compton scattering, especially 
with seed photons from broad line region. In the model fitting to SED including 7-ray emission, the 
external Compton process have been taken into account in combination with SSC (Ghisellini et al. 1998; 
Ghisellini & Tavecchio 2008). The measurements of broad emission lines can be used to constrain the 
energy density of BLR in jets, although it also depends on the location of jet region. On the other 
hand, the measurements of prominent broad emission lines, especially H/3, Mg II and C IV, enable us 
to measure the black hole mass and the bolometric luminosity using various empirical relations, from 
which the accretion mode can be investigated. 

The Gamma ray Large Area Space Telescope [GLAST) was launched on 11 June 2008, and 
renamed the Fermi Gamma Ray Space Telescope shortly after entering its scientific operating mis- 
sion, which began on 1 1 August, 2008. Abdo et al. (2009a) have presented a list of 1 16 bright, > lOcr 
sources at |&| > 10° taken from the list of bright sources (Abdo et al. 2009b) observed with the Fermi 
Large Area Telescope (LAT) in its initial three-month observing period extending from August 4 to 
October 30 of 2008. Of these sources, 106 are associated with blazars with high confidence and com- 
pose the LAT Bright AGN Sample (LBAS). It contains two radio galaxies, namely Centaurus A and 
NGC 1275, and 104 blazars consisting of 57 FSRQs, 42 BL Lac objects, and 5 blazars with uncertain 
classification. The number of low-confidence AGN associations is 1 1 (one source having two possible 
associations - one high and one low confidence). In this work, we search the SDSS photometry and 
spectroscopy data for all these 116 AGNs, and the results focussing on the spectral line analysis are 
shown. The cosmological parameters Hq = 70 km Mpc~^, il„i=0.3, ilA = 0.7 are used throughout 
the paper, and the spectral index a is defined as oc with being the flux density at frequency v. 

2 SDSS DATA ANALYSIS 

The Sloan Digital Sky Survey (SDSSl3 consists of a series of three interlocking imaging and spectro- 
scopic surveys (Legacy, SEGUE, and Supernova), carried out over an eight-year period with a dedicated 
2.5m telescope located at Apache Point Observatory in Southern New Mexico. The seventh data release 
(DR7) from the SDSS represents the completion of this project. The DR7 imaging data cover about 
8423 square degrees of "legacy" sky, with information on roughly 230 million distinct photometric ob- 
jects, and about 3240 square degrees of SEGUE sky, with about 127 million distinct objects (including 
many stars at low latitude). The DR7 spectroscopic data include data from 1802 main survey plates of 
640 spectra each, and cover 8200 square degrees. In addition, DR7 contains 762 "extra" and "special" 
plates. 

Based on the 106 high-confidence association and 11 low-confidence association bright AGN list, 
the photometry data are collected from SDSS DR7 for 28 sources (12 BL Lacs and 16 FSRQs), which 
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is listed in Table [T] Two out of 28 objects are low-confidence association bright AGNs, i.e. BL Lac 
object OFGL J0909.7+0145 identified with SDSS J090939.84+020005.2 {z = 1.575), and FSRQ OFGL 
J1641.4+3939 associated with SDSS J164147.54+393503.3 (z = 0.539). Among these 28 SDSS pho- 
tometric sources, SDSS spectra are available for 20 sources (6 BL Lacs and 14 FSRQs). In below, we 
present the data analysis based on the available SDSS data. 

2.1 Photometry 

SDSS DR7 provides an imaging catalogue of 357 million unique objects, from which the magnitude are 
available at ugriz wavebands with the average wavelengths 355lA, 4686A, 6165A, 748lA and 893lA, 
respectively. As blazars often appear as point-like sources (although the faint host galaxy can be seen for 
BL Lac objects in some case), the PSF magnitude from SDSS are collected for all 28 sources, which is 
determined by fitting a PSF model to the object (see Table 2). The SDSS ugriz magnitudes were firstly 
transferred to AB ugriz magnitudes using the zeropoint offset between the SDSS system and the AB 
system (Oke & Gunn 1983), i.e. u^b = "sdss - 0.04 mag, ^ab ~ ffsDss, '"ab ~ '"sdss, *ab ^ «sdss 
and ZAB = ^SDSS + 0.02 mag. The Galactic extinction are then corrected using the extinction law of 
Schlegel, Finkbeiner & Davis (1998) and the E{B - V) provided in the NASA/IPAC Extragalactic 
Database (NEDfl Finally, flux density can be calculated from magnitudes at each wavebands of ugriz 
as: 

S*^ = 3631 X 10"°* ""Jy (1) 

in which S,y and m are the flux density and Galactic extinction corrected AB magnitude at ugriz wave- 
bands, respectively. The simple power-law 5*,^ = aiy^"" is used to fit the flux density at ugriz wave- 
bands. The spectral index ai, is presented in Table|2l in which the luminosity at 5100A calculated from 
power-law fit is also given. 

2.2 Spectral analysis 

The SDSS spectra cover the wavelength range from 3800 to 9200A with a resolution of about 1800 
- 2200. In a first step, the SDSS spectra for all 20 sources were corrected for the Galactic extinction 
using the reddening map of Schlegel et al. (1998) and then shifted to their rest wavelength, adopting the 
redshift (from NED) listed in Table 1. 

2.2.7 BL Lac objects 

Historically, an object is classified as a BL Lac object based on a flat radio spectrum (ai, < 0.5) and 
an optical spectrum with emission lines weaker than 5 A rest-frame equivalent width, e.g. for the 1-Jy 
radio survey (Stickel et al. 1991; Rector & Stocke 2001). The latter criterion was chosen to separate 
BL Lacs from FSRQ. However, this classification arose controversy. Later on, the absorption feature Ca 
H&K break was proposed in BL Lacs classification in combination with the criteria of no or very weak 
emission lines in the optical spectrum (Marcha et al. 1996; see also Landt, Padovani & Giommi 2002). 
In this work, we follow the BL lacs identification in Abdo et al. (2009a) (see their Tables 1 and 2). 

The SDSS spectra of 6 BL Lacs are featureless, and the Ca H&K break is not evident. Therefore, 
we use a single power-law to fit the spectral continuum, and the contribution of host galaxy is ignored. 
The spectral index is shown in Table 3, as well as the luminosity at 5100A calculated from the fitted 
power-law. 

2.2.2 FSRQs 

The spectra of quasars are characterized by various of broad and narrow emission lines (Vanden Berk 
et al. 2001). The redshift of our 14 FSRQs ranges from 0.539 to 2.189 (see Table 1), with several broad 
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emission lines apparently seen from SDSS spectra, e.g. H/3, Mg II and C IV lines (see Fig. 1). For our 
FSRQs sample, we ignore the host galaxy contribution to the spectrum, since only very little, if any, 
starlight is observed. In order to reliably measure line parameters, we choose those wavelength ranges 
as pseudo-continua, which are not affected by prominent emission lines, and then decompose the spectra 
into the following three components (see also Chen, Gu & Cao 2009): 

1. A power-law continuum to describe the emission from the active nucleus. The 13 line-free spec- 
tral regions were firstly selected from SDSS spectra covering 1 140A to 5630A for our sample, namely, 
1140A-1150A, 1275A-1280A, 1320A-1330A, 1455A-1470A, 1690A- 1700A, 2160A-2180A, 
2225 A - 2250A, 3010A - 3040A, 3240A - 3270A, 3790A - 38 lOA, 4210A - 4230A, 5080A - 5 lOOA, 
5600A - 5630A (Vanden Berk et al. 2001; Forster et al. 2001). Depending on the source redshift, the 
spectrum of individual quasar only covers some of them, from which the initial single power-law are 
obtained for each source. 

2. An Fe II template. The spectra of our sample covers UV and optical regions, therefore, we adopt 
the UV Fe II template from Vestergaard & Wilkes (2001), and optical one from Veron-Cetty et al. 
(2004). For the sources with both UV Fe II and optical Fe II lines prominent in the spectra, we connect 
the UV and optical templates into one template covering the whole spectra (see also Chen et al. 2009). 
In the fitting, we assume that Fe II has the same profile as the relevant broad lines, i.e. the Fe II line 
width usually was fixed to the line width of broad H/3 or Mg II or C IV, which in most cases gave a 
satisfied fit. In some special cases, a free-varying line width was adopted in the fitting to get better fits. 

3. A Balmer continuum generated in the same way as Dietrich et al. (2002) (see also Chen et al. 
2009). Grandi (1982) and Dietrich et al. (2002) proposed that a partially optically thick cloud with a 
uniform temperature could produce the Balmer continuum, which can be expressed as: 

F^^^ = FbeBx{T,){1 - e-^'); (A < Abe) (2) 

where Fbe is a normalized coefficient for the flux at the Balmer edge (Abe = 3646A), Bx{Te) is the 
Planck function at an electron temperature T^, and t\ is the optical depth at A and is expressed as: 

ta = tbe(t^) (3) 
Abe 

where tbe is the optical depth at the Balmer edge. There are two free parameters, Fbe and tbe- 
Following Dietrich et al. (2002), we adopt the electron temperature to be Tg = 15, 000 K. 

The modeling of above three components is performed by minimizing the in the fitting process. 
The final multicomponent fit is then subtracted from the observed spectrum. The fitted power-law, Fe II 
lines, Balmer continuum and the residual spectra for each source are shown in Fig.[T] The Fe II fitting 
windows are selected as the regions with prominent Fe II line emission while no other strong emission 
lines, according to Vestergaard & Wilkes (2001) and Kim et al. (2006). The fitting window around the 
Balmer edge (3625 — 3645A) is used to measure the contribution of Balmer continuum, which extends to 
Mg II line region. The Balmer continuum is not considered when 3625 — 3645A is out of the spectrum. 

The broad emission lines were measured from the continuum subtracted spectra. We mainly focused 
on several prominent emission lines, i.e. H/3, Mg II, C IV. Generally, two gaussian components were 
adopted to fit each of these lines, indicating the broad and narrow line components, respectively. The 
blended narrow lines, e.g. [O III] AA4959, 5007A and [He II] A4686A blending with H/3 were included 
as one gaussian component for each line at the fixed line wavelength. The minimization method was 
used in fits. The line width FWHM, line flux of broad H/3, Mg II and C IV lines were obtained from the 
final fits for our sample, and are listed in Table 4. The spectral fitting for each source are shown in Fig. 

m 

2.3 7\fbh and Lblr 

There are various empirical relations between the radius of broad line region (BLR) and the continuum 
luminosity, which can be used to calculate the black hole mass in combination with the line width 
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FWHM of broad emission lines. However, there are defects when using the continuum luminosity to 
estimate the BLR radius for blazars since the continuum flux of blazars are usually Doppler boosted 
due to the fact that the relativistic jet is oriented close to the line of sight. Alternatively, the broad line 
emission can be a good indicator of thermal emission from accretion process. Therefore for our FSRQs 
sample, we estimate the black hole mass by using the empirical relation based on the luminosity and 
FWHM of broad emission lines. According to the source redshift, we use various relations to estimate 
the black hole mass: Vestergaard & Peterson (2006) for broad H/3 (see also Wu et al. 2004) and Kong et 
al. (2006) for broad Mg II and C IV lines. 

For the sources with available FWHM and luminosity of broad H/3, the method to calculate Afen 
is given by Vestergaard & Peterson (2006): 

In addition, Kong et al. (2006) presented the empirical formula to obtain the black hole mass using 
broad Mg II and C IV for high redshift sources as follows, 

Me„(M. „, . . .„« (^)° Me, 

In the redshift range of our 14 FSRQs, A/bh can be estimated using two of above relations for 7 
sources, in which two AIbh values are consistent with each other within a factor of three. Therefore, we 
adopted the average value of two estimates. 

In this work, we estimate the BLR luminosity Xblr following Celotti et al. (1997) by scaling the 
strong broad emission lines H/3, Mg II and C IV to the quasar template spectrum of Francis et al. (1991), 
in which Lya is used as a reference of 100. By adding the contribution of Ha with a value of 77, the 
total relative BLR flux is 555.77, of which H/3 22, Mg II 34, and C IV 63 (Celotti et al. 1997; Francis et 
al. 1991). From the BLR luminosity, we estimate the bolometric luminosity as Lboi = IOLblr (Netzer 
1990). The black hole mass, BLR luminosity and the Eddington ratio Lboi/iEdd are shown in Table 5. 

3 RESULTS AND DISCUSSION 

In Abdo et al. (2009a), the combination of the figure of merit approach (FoM) and positional association 
methods yields a number of 106 high-confidence (probability P > 0.90) associations (constituting the 
LBAS) and 11 low confidence (0.40 < P < 0.90) associations. These 11 low-confidence association 
Bright AGNs consist of 3 FSRQs, 2 BL Lacs, and 6 sources with uncertain classification, of which all 
have a flat radio spectrum but there are no other data reported in literature allowing to classify them 
either as FSRQs or BL Lacs. For 7-ray source OFGL J1034.0H-6051, two radio associations were found 
by the FoM method, one with very high probability {P = 0.99) and one with lower, but still significant, 
probability (P = 0.79). Although the high-probability source likely dominates the 7-ray emission, it 
is entirely plausible that the low-probability source contributes non-negligibly to the total 7-ray flux 
(Abdo et al. 2009a). In this work, the SDSS data is only available for the high-confidence association 
(see Table 1), from which the emission line measurements are given in Table 4 (also see Table 5). As 
stated in Abdo et al. (2009a), the early results from the first three months of the science mission of 
the Fermi Gamma ray Space Telescope demonstrate its exceptional capabilities to provide important 
new knowledge about 7-ray emission from AGNs and blazars. As the Fermi-LALT data accumulate, 
many more AGNs at lower flux levels will likely be detected helping to improve our understanding of 
supermassive black holes (Abdo et al. 2009a). As predicted from AGN feedback by Wang (2008), there 
is extended faint gamma-ray fuzz around radio-quiet quasars, or red quasars. Fermi Gamma ray Space 
Telescope is anticipated to explore the fuzz more feasibly. There is no doubt that the detection of the 
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fuzz will provide the most powerful diagnostic to the evolutionary chains of galaxies and quasars from 
details of radiative feedback (see Wang 2008). 

Radio-loud quasars, particularly FSRQs, differ from radio quiet AGNs in the contribution of non- 
thermal jet emission in optical continuum emission, in addition to the thermal emission from accretion 
disk. Especially, the jet emission can be dominant due to the beaming effect when jet is moving towards 
us with small viewing angle. Indeed, Liu, Jiang & Gu (2006) found that the 5100A continuum luminosity 
for FSRQs generally exceed that expected from the relationship between the 5100A continuum and 
broad H/3 luminosity for radio quiet AGNs (see also Gu, Chen & Cao 2009). To inspect our FSRQs 
in a similar way, we investigate the relationship between the 3000A continuum and broad Mg II line 
luminosity, since Mg II line is available and measured for most sources. The calibration of this relation is 
firstly performed for a sample of radio quiet AGNs with reverberation mapping black hole mass (Kaspi 
et al. 2000), by collecting the relevant data in Kong et al. (2006). The ordinary least-square bisector 
Unear fit gives: 

log XL^^^^^^ = (0.996 ± 0.048) log Lug ii + (1.895 ± 2.024) (7) 

which is plotted in Fig. 2, as well as the data points of radio quiet AGNs. 

The relationship between the broad Mg II luminosity and luminosity at 3000A for our FSRQs is 
shown in Fig. 2. When Mg II line is not available, we calibrate C IV line to Mg II line adopting the 
relative flux (Mg 11=34, C IV=63) of the relevant lines of the composite spectrum of Francis et al. 
(1991). The luminosity at 3000A is calculated from the power-law fit on spectral continuum. From Fig. 
2, it can be seen that the 3000A luminosity of most sources lies above the XL = — Lme ii relation of 

■' A3000A '= 

radio quiet AGNs. However, the deviation of XL » is within one order of magnitude. The average 

^ A3000A ° ° 

value of Ai^^^^^^/^Mg ii is (log (XL ^^^^^j^/ Lug ii)) = 1-73 ± 0.24 for radio quiet AGNs, while it 
is (log (AL » /Lm. ii)) = 2.21 ± 0.27 for our FSRQs. Therefore, the XL - of our FSRQs is, 

^ A3000A' ^ A3000A ^ 

on average, higher than that of radio quiet AGNs by a factor of three at given LMg ii, which is most 
Ukely due to the non-thermal jet emission. As expected for FSRQs, this implies that the nonthermal 
jet emission in FSRQs generally can be dominant over thermal emission, with the former likely being 
Doppler boosted. Adopting a Doppler factor of 5 = 10 typical for FSRQs (Ghisellini et al. 1998; Jiang, 
Cao & Hong 1998; Gu, Cao & Jiang 2009) and assuming the corresponding Doppler boosting (5^+" 
(corresponding to a continuous jet, a is the spectral index) in jet emission, the intrinsic jet emission in 
optical band, however, can be much smaller than the thermal emission (mainly from accretion disk) for 
our FSRQs. However, the dependence of Doppler factor on the frequency was proposed (Zhang, Fan 
& Cheng 2002), with a lower Doppler factor ~ 5 claimed at optical region. In contrast, Chen et al. 
(2009) shown that the thermal emission from accretion disk and host galaxy can be dominant in 67 out 
of 185 FSRQs (- 36%) selected from SDSS DR 3 quasar catalogue. As Landt et al. (2008) claimed, 
the flat radio spectrum does not prove the jet-dominance in continuum emission. The high radio core- 
dominance, in addition to flat radio spectrum, are proposed to find jet-dominant FSRQs. We note that 
the XL^^^^^j^ — LMg II relation of radio-quiet AGNs is established from relatively low luminosity, which 
does not cover the luminosity range of our FSRQs. The further investigations are needed. 

It can be seen from Table |5] that the black hole masses of 14 FSRQs ranges from IQ^-'^Mq to 
10^-^Mq, with most of sources larger than lO^M©. The BLR luminosity varies from 10^^ '^^erg s^^ 
to lG^^ °'*erg s^^. The Eddington ratio Lboi/LEdd ranges from 10^^-^ to ~ 1. This implies that the 
optically thin, geometric thick accretion disk (Shakura & Sunyaev 1973) may exist in these FSRQs. 
This is consistent with the radio-loud quasars in Wang, Ho & Staubert (2003), of which the accretion 
rates, M w 0.01 — 1 times the Eddington rate, suggesting that most of sources possess standard optically 
thick, geometrically thin accretion disks. In contrast, the accretion mode in BL Lacs is not homogenous, 
with standard thin disks probably in luminous BL Lacs, and the advection-dominated accretion flow in 
less luminous ones (e.g. Cao 2002). Moreover, the evolutionary sequence of blazars, from FSRQs to low- 
energy-peaked BL Lac object to high-energy-peaked BL Lac object, is suggested along the decreasing 
accretion rate (e.g. Cao 2002; Cavaliere & D'Elia 2002; Wang et al. 2003). While we estimate the 
BLR luminosity from multiple emission lines following Celotti et al. (1997), it is argued that the BLR 
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luminosity can be better derived from single line than multiple lines due to the very large scatter in the 
ratio of some key lines (e.g. C IV and Mg II) with respect to Ly a (Wang, Luo & Ho 2004). Following the 
suggestion of Wang et al. (2004), for seven of our FSRQs with multiple lines measurements (see Table 
4), we re-estimated the BLR luminosity using only H/3 line, or Mg II when H/3 is not available, which 
is listed in Table 5. We find that the BLR luminosity from single line is all consistent with that from 
multiple lines within a factor of two, except for SDSS J131028. 66+322043. 7 with the former smaller 
than the latter by a factor of 4.7 (see Table 5). This implies that the relative flux of H/3, Mg II and C IV, 
at least for these seven FSRQs, generally follow the Francis et al. (1991) composite spectrum, although 
they are 7-ray and strong radio sources. Apparently, using single line does not significantly alter the 
Eddington ratio estimates. 

Our BLR measurements can be used in the further investigations relevant to Fermi 7-ray emis- 
sion. For example, the model fit to SEDs including external Compton process and SSC will be present 
for LBAS sources in a future paper (Li et al., in preparation), in which the BLR measurements in this 
paper can be used to constrain the energy density of BLR emission in jets. Moreover, to investigate 
the mechanism of 7-ray emission, the work of Fan et al. (2006), will be re-investigated for LBAS 
sources by using our BLR measurements and collecting the BLR measurements from literatures (Fan et 
al., in preparation). As a matter of fact, the BLR data for five FSRQs (SDSS J082447.24+555242.6, 
SDSS J094857.33+002225.5, SDSS J101603. 14+051302.3, SDSS J155332.69+125651.7 and SDSS 
J222940. 09-083254. 5) are not available from literatures, therefore, are measured and provided in the 
first time in this paper (see Table 4). Despite extensive literature search, the BLR measurements are 
still not available for some LBAS sources, therefore, the spectroscopy observations are needed since the 
external Compton process from BLR seed photons are thought to be primary for 7-ray emission (Fan et 
al. 2006). 
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Table 1 Fermi-detected AGNs in SDSS 
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LAT Name 


Other Name 


SDSS Name 


Class 


z 


SDSS data 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


OFGL J0050.5-0928 


TWT'C^ f\C\ A O C\f\^ 

PKS 0048-097 


J005 04 1.31-092905.1 


BLLac 


0.537 


PS 


OFGL J0738.2+1738 


PKS 0735+178 


J073807. 39+17421 8.9 


BLLac 


0A24 


P 


OFGL J08 18.3+4222 


OJ +425 


J08 18 15. 99+422245. 4 


BLLac 


0.530 


PS 


OFGL J0824.9+5551 


TXS 0820+560 


J082447. 24+555242. 6 


FSRQ 


1.417 


PS 


OFGL J0909.7+0145^ 


PKS 0907+022 


J090939.84+020005.2 


BLLac 


1.575 


PS 


OFGL J092 1.2+4437 


RGB J0920+446 


J092058.46+444154.0 


FSRQ 


2.190 


PS 


OFGL J0948.3+0019 


PMN J0948+0022 


J094857. 33+002225. 5 


FSRQ 


0.585 


PS 


OFGL J0957.6+5522 


4C +55.17 


J095738. 19+552257.7 


FSRQ 


0.896 


PS 


OFGL J1015.2+4927 


lES 1011+496 


J101504. 14+492600.6 


BLLac 


0.212 


PS 


OFGL J1015.9+0515 


PMN J1016+0512 


J101603. 14+05 1302.3 


FSRQ 


1.715 


PS 


OFGL J1034.0+6051 


S4 1030+61 


J103351.42+605107.3 


FSRQ 


1.401 


PS 


OFGL J1057.8+0138 


PKS 1055+018 


J105829.60+013358.8 


FSRQ 


0.888 


PS 


OFGL J1058. 9+5629 


RX J10586+5628 


J 105837.73+5628 11.1 


BLLac 


0.143 


PS 


OFGL Jl 104.5+3811 


Mrk421 


J 110427.31+381231.7 


BLLac 


0.030 


p 


OFGL Jl 159.2+2912 


4C +29.45 


Jl 15931.84+291443.8 


FSRQ 


0.724 


PS 


OFGL J1218.0+3006 


B2 1215+30 


J121752.08+300700.6 


BLLac 


0.130 


p 


OFGL J1221.7+2814 


WCom 


J122131.69+281358.4 


BLLac 


0.102 


PS 


OFGL J 1229. 1+0202 


3C 273 


J122906.69+020308.5 


FSRQ 


0.158 


p 


OFGL J13 10.6+3220 


B2 1308+32 


J131028.66+322043.7 


FSRQ 


0.997 


PS 


OFGL J1427. 1+2347 


PKS 1424+240 


J142700.39+234800.0 


BLLac 




p 


OFGL J1504.4+1030 


PKS 1502+106 


J150424.98+102939.1 


FSRQ 


1.839 


PS 


OFGL J 1522.2+3 143 


TXS 1520+319 


J152209.99+314414.3 


FSRQ 


1.487 


p 


OFGL J1553.4+1255 


PKS 1551 + 130 


J155332.69+125651.7 


FSRQ 


1.290 


PS 


OFGL J1555. 8+1 110 


PG 1553+11 


J155543.04+1 11 124.3 


BLLac 


0.360 


p 


OFGL J1635. 2+3809 


4C +38.41 


J163515.50+380804.4 


FSRQ 


1.813 


PS 


OFGL J 164 1.4+3939^ 


B3 1640+396 


J164147.54+393503.3 


FSRQ 


0.539 


PS 


OFGL J1653.9+3946 


Mrk501 


J165352.21+394536.6 


BLLac 


0.033 


p 


OFGL J2229.8-0829 


PHL 5225 


J222940.09-083254.5 


FSRQ 


1.560 


PS 



Notes: Col. (1): source LAT name. ^ represents low-confidence association AGNs. Col. (2): source other 
name. Col. (3): source SDSS name. Col. (4): class. Col. (5): redshift from NED. Col. (6): available SDSS data 
— P for photometiy and S for spectroscopy. 
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Table 2 SDSS photometry 



SDSS Name u g r i z E(B - V) a„ log L 

(erg s ) 



(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


JUUjU41, 


.J L-yjyZyVj.l 


16 


.77 


16 


.39 


16 


.05 


1 ^ on 
1 j.SU 


1 ^ ^7 

1 J.J / 


n mo 


1 


.03 


4j.yj 


jV / joV / . 


.jy+l /4Z10.7 


16 


.45 


15, 


.91 


15 


.49 


1 J. 1 / 


1/1 QQ 
14. Ko 


U.Ujj 


1 


.38 


4j.yj 


TAO 1 O 1 C 


QQ_i_/1 000/1 ^ A 


19, 


.20 


18, 


.61 


18, 


.08 


17 71 

1 / . / 1 


1 7 

1 /.JO 


U.UOj 


1 


.55 


/I ^ 1 C 
4j. lo 


JUoZ44 / , 


0/1_l_'^^'^0/10 A 


18, 


.52 


18, 


.36 


18, 


.08 


1 7 QO 


1 7 QC 


U.UOj 


0, 


.43 


A ^ QO 


J090939, 


.84+020005.2^ 


20, 


.37 


19, 


.88 


19, 


.53 


19.18 


18.98 


0.031 


1 


.24 


45.76 


J092058, 


.46+444154.0 


18, 


.68 


18, 


.06 


17, 


.91 


17.72 


17.42 


0.021 


1 


.02 


46.68 


J094857, 


.33+002225.5 


18, 


.96 


18, 


.60 


18, 


.44 


18.19 


18.18 


0.079 


0, 


.48 


45.05 


J095738, 


.19+552257.7 


18, 


.26 


17, 


.92 


17, 


.57 


17.38 


17.20 


0.009 


0, 


.99 


45.81 


T 1 n 1 c n A 

J101504, 


1 A . A C\'~\ /' r\f\ f 

.14+492600.6 


15 


.71 


15 


.43 


15 


.23 


15.1 1 


14.96 


0.012 


0, 


.64 


45.28 


Ti ni 


1 /I-lH^ 1 ^no ^ 


20, 


.50 


19 


.88 


19, 


.48 


1 Q n7 


1 Q 7'^ 
18. / J 


n noA 
u.uzo 


1 


.58 


/I A nn 
40. UU 


Ti n^^^ 1 




19 


.31 


18, 


.88 


18, 


.46 


1 Q 1 7 
lo. 1 / 


1 7 

1 /.Vj 


n ni n 
u.uiu 


1 


.30 


/I A i\A 
40. U4 


TI n^SOQ 


.0U+U13 J jo.o 


18, 


.84 


18, 


.30 


17, 


.85 


1 7 Al 


17 10 
1 / . IZ 


n no7 
u.uz / 


1 


.56 


4J.OJ 


J105837, 


.73+562811.1 


16, 


.92 


16, 


.50 


16, 


.00 


15.73 


15.49 


0.012 


1 


.36 


44.57 


J110427, 


.31+381231.7 


14, 


.80 


13, 


.81 


13, 


.10 


12.81 


12.54 


0.015 


2, 


.14 


44.19 


J115931 


.84+291443.8 


18, 


.75 


18, 


.21 


18, 


.08 


17.95 


17.68 


0.019 


0, 


.85 


45.39 


J121752, 


.08+300700.6 


15, 


.92 


15, 


.49 


15, 


.14 


14.92 


14.76 


0.024 


1 


.03 


44.85 


J122131 


.69+281358.4 


15, 


.84 


15, 


.39 


15, 


.05 


14.78 


14.58 


0.023 


1 


.13 


44.67 


J 122906, 


.69+020308.5 


13, 


.87 


12, 


.99 


12, 


.88 


12.64 


13.24 


0.021 


0, 


.61 


45.88 


J131028, 


.66+322043.7 


18, 


.40 


18, 


.02 


17, 


.61 


17.34 


17.12 


0.014 


1 


.20 


45.97 


J142700, 


.39+234800.0 


15, 


.25 


14, 


.87 


14, 


.55 


14.31 


14.12 


0.059 


0, 


.88 




J 150424, 


.98+102939.1 


18, 


.90 


18, 


.52 


18, 


.22 


17.84 


17.59 


0.032 


1 


.15 


46.45 


J 152209, 


.99+314414.3 


20, 


.62 


20, 


.18 


19 


.79 


19.51 


19.43 


0.024 


1 


.10 


45.52 


J155332, 


.69+125651.7 


17, 


.50 


17, 


.40 


17, 


.24 


17.20 


17.17 


0.042 


0, 


.16 


46.14 


J155543, 


.04+111124.3 


15, 


.08 


14, 


.74 


15, 


.40 


14.13 


13.94 


0.052 


0, 


.84 


46.10 


J163515, 


.50+380804.4 


17, 


.77 


17, 


.67 


17, 


.63 


17.38 


17.26 


0.011 


0, 


.42 


46.40 


J164147, 


.54+393503.3^ 


20, 


.09 


19, 


.65 


19, 


.64 


19.23 


19.16 


0.014 


0, 


.81 


44.50 


J165352, 


.21+394536.6 


15, 


.36 


13, 


.85 


13, 


.04 


12.62 


12.32 


0.019 


2, 


.86 


44.24 


J222940, 


.09-083254.5 


17, 


.87 


17, 


.44 


17, 


.07 


16.74 


16.48 


0.051 


1 


.16 


46.74 



Notes: Col. (1): SDSS name. ^ represents low-confidence association AGNs. Cols. (2) - (6): SDSS five-band 
PSF magnitude ugriz. Col. (7): E{B — V) from NED. Col. (8): the spectral index from power-law fit on 
u, g, r, i, z. Col. (9): luminosity at 5100A. 



Table 3 BL Lacs with SDSS spectra 



SDSS Name ax log L 

^ 5 100 A 

(erg s ) 



(1) 


(2) 


(3) 


J005041.31-092905.1 


0.69 


46.19 


J081815.99+422245.4 


0.79 


45.00 


J090939.84+020005.2^ 


1.05 


46.63 


J101504.14+492600.6 


1.32 


45.37 


J105837.73+562811.1 


0.87 


44.69 


J122131.69+281358.4 


0.86 


44.85 



Notes: Col. (1): SDSS name. ^ represents low-confidence association AGNs. Col. (2): the spectral index ax 
from power-law fit to SDSS spectra. Col. (3): luminosity at 5100A. 
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Table 4 Broad emission line measurements for FSRQs 
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C TV 

V_ i V 














" 3000A 




FWHM 


log F 


FWHM 


log F 


FWHM 


log F 






(kms ) 




(kms ) 




(kms ) 






(erg s ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


J082447.24+555242.6* 






5119 


-13.27 






1.39 


46.88 


J092058.46+444154.0 


8788 


-12.68 










1.18 


47.69 


J094857.33+002225.5* 






4743 


-13.89 


2209 


-14.15 


1.28 


45.50 


J095738. 19+552257.7 






6709 


-13.70 






0.93 


46.30 


J101603. 14+05 1302.3* 






5268 


-13.66 






0.75 


46.99 


J103351.42+605107.3 






7607 


-14.04 






0.43 


46.34 


J105829.60+013358.8 






6214 


-13.48 






0.61 


46.38 


Jl 15931.84+291443.8 






5008 


-13.24 


4636 


-13.56 


0.52 


46.67 


J131028.66+322043.7 






5278 


-13.36 


3633 


-14.40 


0.76 


46.51 


J150424.98+102939.1 


6627 


-12.90 


5172 


-13.35 






1.69 


47.44 


J155332.69+125651.7* 






3999 


-13.34 






1.28 


47.24 


J1635 15.50+380804.4 


7452 


-12.52 


5144 


-13.04 






1.08 


47.52 


1164147.54+393503.3^ 






7064 


-13.44 


7731 


-13.72 


1.84 


45.11 


J222940.09-083254.5* 


7753 


-12.69 


3689 


-13.27 






1.53 


46.81 



Notes: Col. (1): SDSS name. * for sources with BLR data not available in literatures and newly measured 
in tliis paper ^ represents low-confidence association AGNs. Col. (2): FWHM of broad C IV. Col. (3): flux of 
broad C IV in units of org cm^^. Col. (4): FWHM of broad Mg II. Col. (5): flux of broad Mg II in units 
of erg cm~^. Col. (6): FWHM of broad H/3. Col. (7): flux of broad H/3 in units of erg s^^ cm^^. Col. 
(8): the continuum spectral index ax from SDSS spectra. Col. (9): luminosity at 3000A. 



Table 5 Black hole mass and BLR luminosity 


SDSS Name 


log i\/BH 


log I/BLR 


log Lboi 


log -£/bol/ L-Edd 


log l'blh 




(Mq) 


(erg s"^) 


(ergs^^) 




(ergs^^) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


J082447.24+555242.6 


9.49 


46.04 


47.04 


-0.56 




J092058.46+444 154.0 


9.88 


47.04 


48.04 


0.05 




J094857.33+002225.5 


8.26 


44.44 


45.44 


-0.93 


44.39" 


J095738. 19+552257.7 


9.20 


45.11 


46.11 


-1.20 




J101603. 14+051302.3 


9.78 


46.00 


47.00 


-0.89 




J103351.42+605107.3 


9.39 


45.25 


46.25 


-1.25 




J105829.60+013358.8 


9.25 


45.32 


46.32 


-1.04 




Jl 15931.84+291443.8 


9.11 


45.30 


46.30 


-0.92 


45.22" 


J131028.66+322043.7 


8.94 


45.39 


46.39 


-0.66 


44.72" 


J150424.98+102939.1 


9.50 


46.36 


47.36 


-0.25 


46.23'' 


J155332.69+125651.7 


9.19 


45.88 


46.88 


-0.42 




J1635 15.50+380804.4 


9.74 


46.71 


47.71 


-0.14 


46.53'' 


J164147.54+393503.3^ 


9.19 


44.80 


45.80 


-1.50 


44.74" 


J222940.09-083254.5 


9.40 


46.36 


47.36 


-0.15 


46.13'' 



Notes: Col. (1): SDSS name. ^ represents low-confidence association AGNs. Col. (2): the black hole mass. 
Col. (3): BLR luminosity. Col. (4): bolometric luminosity. Col. (5): the Eddington ratio. Col. (6): BLR luminos- 
ity estimated from single line for H/3, ^ for Mg II). 
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400 - 





SDSS J092058. 46 + 444154,0 
z = 2.189 




2400 2600 2800 3000 3200 3400 3600 
Rest Wavelength (A) 



T 200 
°< 

'a 1 50 

o 

'cfi 1 00 

bC 
u 

i> oO 
o 

-50 



SDSS J094857. 33 + 002225, 5 

z = 0.585 2 




o 200 



100 r 



1500 2000 2500 

Rest Wavelength (A) 

SDSS J09573'8. 19 + 552257.7 
z = 0.896 




2500 5000 3500 4000 4500 5000 
Rest Wavelength (A) 



SDSS J101603. 14 + 051302. 3 
z=1 .71 5 



2500 3000 3500 4000 4500 

Rest Wavelength (A) 



500 - 



I 200 fW»i;.,.,,v,.H,;^ 




SDSS J103351. 42 + 605107. 3 
z=1 .401 




2000 2500 
Rest Wavelength (A) 



3000 2400 2600 2800 3000 3200 3400 3600 3800 
Rest Wavelength (A) 



1500 - 



£r 500 



SDSS J105829. 60 + 013358. 8 
feu, , z = 0.888 



1 000 
-< 800 



600 
'» 400 
; 200 

- 

-200 



SDSS J1 15931.84 + 291443.8 
z = 0.724 



2500 3000 3500 4000 4500 

Rest Wavelength (A) 



2500 3000 3500 4000 4500 5000 
Rest Wavelength (A) 



Fig. 1 The spectral analysis. The source SDSS name and redshift are presented in each panel 
in sequences of increasing source R.A.. For each source panel, the top and middle black lines 
are the original and the continuum-subtracted spectrum, respectively. The top dashed red line 
represents the power-law continuum. The middle blue lines are the individual line components 
in multi-line spectral fitting, and the middle solid red line is integrated lines fitting. The green 
spectral region in the continuum-subtracted spectrum is used to fit Fe II emission and power- 
law continuum, and the bottom green line is Fe II emission, which is shifted downwards with 
arbitrary unit for the sake of presentation. When necessary, the Balmer continuum is indicated 
as solid blue lines in the bottom of each panel, which begins at the Balmer edge 3646A. 
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Fig. 1 continue. 
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40 I 1 1 1 1 1 1 1 

39 40 41 42 43 44 45 46 

Fig. 2 Luminosity at 3000A versus broad Mg II lu- 
minosity. The open circles are our FSRQs, while open 
triangles represent the radio quiet AGNs in Kong et al. 
(2006). The solid line is the OLS bisector linear fit to 
radio quiet AGNs in Kong et al. (2006), XL <■ = 

^ & \ A3000A 

78.5 i^M^^n (^^^ equation |7ll. The dashed line repre- 
sents the deviation from the solid line by one order 
of magnitude in XL % which is used to indicate 

° A3000A 

the deviation of our FSRQs from the relation for radio 
quiet AGNs. 



